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As a result of an extensive observational campaign targeting the Virgo cluster, we obtained 
integrated (drift-scan mode) optical spectra and multiwavelength (UV,U,B,V,H) photometry 
for 124 and 330 galaxies, respectively, spanning the whole Hubble sequence, and with m p < 
16 (M p < —15). These data were combined to obtain galaxy Spectral Energy Distributions 
(SEDs) extending from 2000 to 22000 A. By fitting these SEDs with synthetic ones derived using 
Bruzual & Chariot population synthesis models we try to constrain observationally the Star 
Formation History (SFH) of galaxies in the rich cluster of galaxies nearest to us. Assuming a 
Salpeter IMF and an analytical form for the SFH, the fit free parameters arc: the age (T) of the 
star formation event, its characteristic time-scale (r) and the initial metallicity (Z). In this work 
we test the (simplistic) case where all galaxies have a common age T=13 Gyr, exploring a SFH 
with "delayed" exponential form (which we call "a la Sandage"), thus allowing for an increasing 
SFR with time. This SFH is consistent with the full range of observed SEDs, provided that the 
characteristic time-scale r is let free to vary between 0.1 (quasi instantaneous burst) and 25 Gyr 
(increasing SFR) and Z between 1/50 and 2.5 Z©. Elliptical galaxies (including dEs) are best 
fitted with short time scales (r ~ 3 Gyr) and metallicity varying between 1/5 and 1 Z©. The 
model metallicity is found to increase as a function of H band luminosity. Spiral galaxies require 
that both r and metallicity correlate with H band luminosity: low mass Im+BCD have sub-solar 
Z and r > 10 Gyr, whereas giant spirals have solar metallicities and r ~ 3 Gyr, consistent with 
elliptical galaxies. Moreover we find that the SFH of spiral galaxies in the Virgo cluster depends 
upon the presence at their interior of fresh gas capable of sustaining the star formation. In fact 
the residuals of the r vs. Ljj relation depend significantly on the HI content. HI deficient galaxies 
have shorter (up to a factor of 4) r (truncated SFH) than spirals with normal HI content. 
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Introduction 



1 Based on observations collected at the Observatoire de 
Haute Provence (OHP) (France), operated by the CNRS, 
France, and at the European Southern Observatory (Chile) 
(programme 66.B-0026) 



One of the hottest yet unsettled issues of obser- 
vational cosmology is the tracing of galaxy evolu- 
tion over a large interval of look-back-time, in or- 



der to shed light on the processes that led to their 
formation. Currently debated formation models 
can be broadly divided into two categories: the 
hierarchical (Cole et al. 1994) or the monolithic 
collapse (Eggen, Lynden-Bell & Sandage 1962; 
Sandage 1986) models (see a review by Elmegreen 
2001). 

This issue can only now be pursued observation- 
ally, with the present generation of 10 m tele- 
scopes. Photometric and spectro-photometric ob- 
servations of selected samples of galaxies spanning 
a large redshift interval can provide us with the 
time dependence of their structural parameters, 
their dynamics, their gaseous content and metal 
enrichment (Kauffmann & Chariot 1998) as well 
as their star formation histories. 
An alternative approach to the issue (eg. Bell & 
de Jong 2000) is based on the assumption that 
galaxies at z—0 retain some memory of their past, 
which can be unveiled by models (e.g. Bruzual & 
Chariot, 1993 population synthesis models). Ob- 
servations at z=0 would provide the boundary 
conditions for the modeling process. 
The fundamental ingredient of both approaches 
is a robust observational determination of the 
"end-point" of galaxy evolution: i.e. the present 
stage galaxy properties. Even this relatively sim- 
ple point is far from settled, since the properties 
of local galaxies are not yet sufficiently known and 
understood. For example, it has been shown that 
the present star formation activity (per unit mass) 
increases along the Hubble sequence (Robertsfe 
Haynes 1994; Kennicutt 1998), but very little at- 
tention has been given to the role of total stel- 
lar mass in determining the present evolution of 
galaxies, an issue that we addressed in Gavazzi et 
al. (1996a), Gavazzi & Scodeggio (1996), Boselli 
et al. (2001), and we reiterate in the present work 
(see also Bell & de Jong 2000 who stressed the 
role of the central surface brightness in regulating 
the evolution). 

The properties of the underlying stellar popula- 
tion (and the physical conditions of the interstellar 
medium) can be best addressed using spectro- 
photometric measurements which are unfortu- 
nately only sparsely available. After the pioneer- 
ing work of Kennicutt (1992), marginal effort was 
devoted at deriving the spectral parameters repre- 
sentative of normal galaxies in the local Universe. 
The work of Jansen et al. (2000), who reports 



on ~ 200 spectra of local isolated galaxies, repre- 
sents a significant exception. In the near future 
our knowledge of evolved galaxies will certainly 
benefit from the SIRTF legacy project SINGS 
(http://ircamera.as.arizona.edu/legacy/) by Ken- 
nicutt et al. specifically aimed at obtaining the 
detailed phenomenology of local galaxies through 
all possible observational windows. 
Another aspect that lately has not been given 
enough consideration when studying galaxy for- 
mation and evolution is the role of environment 
in regulating those processes. Clusters of galax- 
ies are " laboratories" where galaxy evolution took 
place in significantly different environmental con- 
ditions from the field, the most evident signature 
of this being the morphology-density correlation 
(Dressier 1980), but other more subtle effects are 
known and others wait to be quantified. 
With the aim of constructing a representative de- 
scription of galaxies in a nearby cluster, we under- 
took a multi-frequency photometric survey of op- 
tically selected galaxies in the Virgo cluster, span- 
ning the broadest possible range in morphological 
type (Ellipticals, spirals, dE, Im and BCD) and 
luminosity (-22 < Mb < -15). Due to its nearness, 
the Virgo cluster offers a unique opportunity for 
carrying out such an analysis, because dwarf (Mb 
< -15) galaxies in this cluster are within reach of 
middle-class telescopes and, owing to the monu- 
mental photographic work of Sandage, Binggeli & 
Tammann which gave origin to the Virgo Cluster 
Catalogue (Binggeli et al. 1985; VCC), for galax- 
ies in this cluster we can rely on a particularly 
accurate morphological classification. 
Our own Ha, optical, ncar-IR, and millimctric 
observations, together with UV (2000 A) and cen- 
timetric data taken from the literature, were col- 
lected in a multifrcquency database. 
Most recently we undertook a spectro-photometric 
project aimed at obtaining an homogeneous spec- 
troscopic data-set for these galaxies. Integrated 
spectra of 125 Virgo galaxies have been obtained 
so far, covering the spectral range from 3600 to 
7000 A, with a resolution of i?=500-1000. The 
spectra will be published in a forthcoming paper 
(Paper II of this series). 

In the present paper we combine the newly ob- 
tained spectro-photometric data with the photo- 
metric observations to map each galaxy spectral 
energy distribution (SED) over a very broad wave- 
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length range. These SEDs are then used with 
Bruzual & Chariot (1993, hereafter B&C) popu- 
lation synthesis models to constrain the star for- 
mation histories of galaxies in the Virgo cluster. 
The sample upon which the present work is based 
is discussed in Section 2. The observations are 
briefly described in Section 3. The construction of 
the galaxy SEDs, by combining photometry with 
spectroscopy, and their corrections are discussed 
in Section 4. The method adopted to fit B&C 
models to the data is described in Section 5. The 
results of this work are given in Section 6 and 
discussed in Section 7. 

2. The sample 

The Virgo Cluster Catalogue (VCC, Binggeli 
et al. 1985) lists 598 bona-fide Virgo cluster mem- 
bers brighter than m p = 16.0. Of these, 552 are 
members spectroscopically confirmed by Binggeli 
ct al. (1993) or by Gavazzi et al. (2000c), while 
the remaining 46 are classified as " possible cluster 
members". Among the 598 galaxies, 312 are of 
late- type (from Sa to Im-BCD). 
The availability of photometric data for these VCC 
galaxies is detailed in Table 1. 
For the purpose of this paper we selected, among 
Virgo members with available photometry, a sub- 
sample of 330 objects (only 4 arc not spectro- 
scopically confirmed members) according to the 
following criteria: 

1) They have the V magnitude plus at least one 
NIR magnitude (J or H or K; H in most cases) and 
at least one blue/UV magnitude (B or U or UV 
or Ha) measurement available. Fig.l provides a 
graphical representation of this selection criterion. 

2) To avoid large errors on total magnitudes that 
would derive from the extrapolation of poorly de- 
fined growth curves, we require B, V and NIR 
magnitudes to be obtained from at least 3 aper- 
ture measurements (only 5 galaxies were rejected 
because of this requirement. This same require- 
ment does not apply to UV magnitudes because 
these are already total magnitudes). 

Among these 330 objects, 124 galaxies with avail- 
able spectra form the " spectro-photometric" sam- 
ple. Galaxies meeting criteria 1) and 2), but 
without spectroscopic observations, constitute our 
" photometric sample" . 

Details of the sample completeness in three bins 
of photographic magnitude and for the individual 




U(J.H,K)n(U(H„.UV,U,B))nv 

Fig. 1. — 330 objects meet the condition for inclu- 
sion in the "photometric sample" (with m p < 16). 
Of these 124, observed spectroscopically, consti- 
tute the " spectro-photometric sample" . 




12 40 12 20 

R.A. (1950) 



Fig. 2. — Sky distribution of the 598 galax- 
ies brighter than m p < 16.0 in the VCC. The 
filled symbols represent galaxies in the spectro- 
photometric sample, the empty ones galaxies in 
the photometric sample, while the tiny symbols 
refer to unobserved galaxies. Elliptical galax- 
ies are given with squares, spiral galaxies with 
Dcfni < 0.5 are triangles, spiral galaxies with 
Defni > 0.5 are circles. Circles are drawn at 2, 4, 
6 deg. projected radial distance from M87 (cross). 
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Table 1 

The completeness in the individual photometric bands for the 598 Virgo members (or for 

the 312 spirals) listed in the vcc with m p < 16 



N (%) 



UV 134/598 (22) 

U 302/598 (51) 

B 393/598 (66) 

V 383/598 (64) 

J 193/598 (32) 

H 355/598 (59) 

K 255/598 (43) 



(spirals) Ha 228/312 (73) 
(spirals) HI 293/312 (94) 



Table 2 

The available photometric measurements in the individual bands of the 
spectro-photometric and photometric samples, in three bins of photographic magnitude 



Band 




N Phot % 


N Phot+Spec % 


N S 


% 


UV 


<16 


57/598 (10) 


63/598 (11) 


120/598 


(21) 




<15 


54/414 (13) 


62/414 (15) 


116/414 


(28) 




<14 


43/248 (17) 


50/248 (20) 


93/248 


(37) 


u 


<16 


161/598 (27) 


105/598 (18) 


266/598 


(45) 




<15 


150/414 (36) 


94/414 (23) 


244/414 


(59) 




<14 


125/248 (50) 


77/248 (31) 


202/248 


(81) 


B or V 


<16 


206/598 (34) 


124/598 (21) 


330/598 


(55) 




<15 


178/414 (43) 


109/414 (26) 


287/414 


(69) 




<14 


136/248 (55) 


86/248 (35) 


222/248 


(90) 


J 


<16 


118/598 (20) 


64/598 (11) 


182/598 


(31) 




<15 


105/414 (25) 


57/414 (14) 


162/414 


(39) 




<14 


89/248 (36) 


51/248 (21) 


140/248 


(57) 


H 


<16 


191/598 (32) 


100/598 (17) 


291/598 


(49) 




<15 


171/414 (41) 


91/414 (22) 


262/414 


(63) 




<14 


134/248 (54) 


77/248 (31) 


211/248 


(85) 


K 


<16 


136/598 (23) 


101/598 (17) 


237/598 


(40) 




<15 


117/414 (28) 


87/414 (21) 


204/414 


(49) 




<14 


96/248 (39) 


72/248 (29) 


168/248 


(68) 


spir. HI 


<16 


114/312 (37) 


86/312 (28) 


200/312 


(65) 




<15 


98/234 (42) 


77/234 (33) 


175/234 


(75) 




<14 


78/150 (52) 


60/150 (40) 


138/150 


(92) 


spir. Ha 


<16 


92/312 (29) 


87/312 (28) 


179/312 


(57) 




<15 


78/234 (33) 


78/234 (33) 


156/234 


(66) 




<14 


62/150 (41) 


60/150 (40) 


122/150 


(81) 
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Table 3 

The completeness of the spectro-photometric and photometric samples in three bins of 

photographic magnitude. 



mag. limit 


N VCC 


with z 


N Phot 


(%) 


N Phot+Spec 


(%) 


N S 


(%) 


(incomplete) 






220 




125 




345 




<16 


598 


552 


206 


(34) 


124 


(21) 


330 


(55) 


<15 


414 


411 


178 


(43) 


109 


(26) 


287 


(69) 


<14 


248 


248 


136 


(55) 


86 


(35) 


222 


(90) 


6<2orA<9<6 fromM87 


















(incomplete) 






85 




108 




193 




<16 


270 


244 


73 


(27) 


108 


(40) 


181 


(67) 


<15 


185 


184 


53 


(29) 


95 


(51) 


148 


(80) 


<14 


107 


107 


28 


(26) 


75 


(70) 


103 


(96) 



Table 4 

The new spectra in the Virgo Sample. 



Obs. run 


Tel. 


Istr. Resolution 


Slit Width 


N. Obs. 








arcsec 




OHP 1998 


1.93m 


Carelec 500 


2.5 


10 


OHP 1999 


1.93m 


Carelec 1000 


2.5 


17 


OHP 2000 


1.93m 


Carelec 1000 


2.5 


26 


OHP 2001 


1.93m 


Carelec 1000 


2.5 


24 


ESO 2001 


3.6m 


EFOSC2 500 


1.5 


39 


WHT 2001 


4.2m 


ISIS 2000 


2.5 


7 






Table 5 








The linear regression analysis 




sample 




regression 


R 


see Fig. 


Ellipticals 


Logr 


-0.086 x LogL H + 1.351 


-0.537 


11 


Ellipticals 


1! II 


= 0.302 x LogL H + 0.599 


0.711 


15 


Spirals 


Logr = 


-0.149 x LogLH + 2.221 


-0.604 


13 


Spirals 


ALogr 


= -0.132 x Derm + 0.071 


-0.529 


16 



5 



photometric bands are given separately for the 
"photometric" and " spectro-photometric" sam- 
ples, and for their sum in Table 2, while the final 
samples and their completeness are illustrated in 
Table 3. Their sky distribution is shown in Fig. 

2. While the global completeness of the "spectro- 
photometric sample" is still rather poor at m p < 
16.0 (20%), it improves significantly (to 40%) over 
the restricted area composed by the annuli with 
angular distance 6 < 2 deg and 4 < 8 < 6 from 
M87 where we concentrated most of our spectro- 
scopic effort (108 spectra). 

In addition to the 330 galaxies selected with the 
above criteria we include in our analysis 15 ob- 
jects with m p > 16.0 (mostly dwarf ellipticals and 
BCDs), out of which 14 have photometry and one 
(VCC636) has a spectrum available. These are: 
VCC 328, 636, 793, 810, 872, 882, 916, 1065, 1173, 
1313, 1352, 1353, 1377, 1420, 1577. This addition 
brings to 345 and 125 the number of galaxies 
in the "photometric" and "spectro-photometric" 
samples, respectively (as listed in the first line of 
Table 3). 

Following Gavazzi et al. (1999b) we assume a dis- 
tance of 17 Mpc for the members (and possible 
members) of Virgo cluster A, 22 Mpc for Virgo 
cluster B, 32 Mpc for objects in the M and W 
clouds, adopting H a = 75 km s _1 Mpc _1 . 

3. The Data 

3.1. Spectro-photometry 

Spectro-photometric measurements of galaxies 
in the Virgo cluster were taken during several runs 
from 1998 to 2001, using various telescopes. The 
detailed presentation of the spectra, including the 
line analysis, will be addressed in a forthcoming 
paper (Gavazzi et al. Paper II, in preparation). 
Here we give some concise information focused 
on the continuum spectral properties. Out of the 
125 spectra used in this work, 2 were taken from 
the catalogue of Kennicutt (1992). Spectra of 7 
BCD galaxies, obtained with the William Hcrschel 
Telescope with a resolution of 2000, were kindly 
provided to us by J.M. Vilchez. All the remaining 
spectra were obtained by us: 77 using the OHP 
1.93m and 39 using the ESO 3.6m telescope (see 
Table 4). 

The observations were taken in " drift mode" : i. e. 



with the slit, generally parallel to the galaxy ma- 
jor axis, drifting over most of the optical surface of 
the galaxy (as in Kennicutt 1992). Spectra taken 
in this way are representative of the entire galaxy, 
and not only of its central regions. All spectra 
cover the wavelength range 3600-7000 A with a 
resolution of 500 (ESO) and 1000 (OHP). Spectra 
were flux-calibrated using several spectrophoto- 
metric standards. However, because of the use of 
drift-scan mode, their absolute flux calibration is 
meaningless, and therefore all spectra were nor- 
malized to their value at A = 5500 A. 

3.2. Photometry 

Total UV, optical (U,B,V) and NIR (J,H,K) 
magnitudes are used to complement the spec- 
troscopic observations in the construction of the 
galaxies SEDs on which this work is based. 
U, B and V photometry is generally derived from 
our own CCD measurements. When these are not 
available it is derived from aperture photometry 
taken from the Longo et al. (1983) catalogue and 
the LEDA database (Prugnel & Heauderau 1998). 
We used only galaxies with at least three indepen- 
dent aperture measurements, after checking for 
foreground star contamination. We compute to- 
tal magnitudes (U,B,V)25, i.e. magnitudes com- 
puted at the optical radius (as in Gavazzi & Boselli 
1996), which are, on average, 0.1 mag fainter than 
the total asymptotic magnitudes. 
NIR data from Nicmos3 observations, obtained 
by us, have been presented in Gavazzi et al. 
(1996b, c), Gavazzi et al. (2000a), Gavazzi et al. 
(2001), Boselli et al. (1997), Boselli et al. (2000). 
In most cases we measured the H band magnitude, 
in many cases the K band one, in a few cases the 
J band one. From these data we derive total mag- 
nitudes H25, determined as described in Gavazzi 
& Boselli (1996). When H25 is not available we 
derive it from K25 using H — K=0.25 mag. H25 
values are converted into total luminosities using: 
logL H = 11.36 - 0.4H 25 + 2LogD (in solar units), 
where D is the distance to the source (in Mpc). 
The assumed photometrical uncertainties are 15 
% for B, V and H and 20 % for U, J and K (see 
also appendix C). 

For 120 of the 345 galaxies in the photometric 
sample (63 in the spectro-photometric sample) UV 
data are available from three balloon experiments: 
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the FOCA (Milliard et al. 1991), FAUST (Lamp- 
ton ct al. 1993) and SCAP (Donas et al. 1987). 
SCAP and FOCA use 150 A wide filters cen- 
tered at 2000 A, while FAUST is centered at 1650 
A (AA = 2501). The published FOCA magni- 
tudes (Donas ct al. 1991, 1995), were recently 
corrected (Donas, private communication) to ac- 
count for better cross-calibration between the var- 
ious experiments. These are total (not aperture) 
UV magnitudes extracted from the photographic 
plates. The FAUST data were converted by us to 
2000 A (using 0.2 mag correction on average, as 
concluded by Deharveng et al. 1994). The quoted 
error on the UV magnitude is 0.3 mag in general, 
but it ranges from 0.2 mag for bright galaxies to 
0.5 mag for weak sources observed in frames with 
larger than average calibration uncertainties. We 
assume a conservative 0.4 mag uncertainty for all 
galaxies to account for all sources of systematic 
errors. 

3.3. Other data 

Ha+[NII] fluxes enter indirectly in the SEDs 
determination, as they provide an estimate of the 
ionizing flux below 912 A (see Section 4.1 and 
Appendix A) . They were measured in our spectra, 
when available, or derived from imaging or aper- 
ture photometry by Kennicutt & Kent (1983), Al- 
moznino & Brosch (1998), Heller et al. (1999), 
Koopmann et al. (2001) and references therein. 
Additional imaging observations of 150 galaxies 
have been recently obtained by us during sev- 
eral runs at the Observatoire d'Haute Provence 
(France), at Calar Alto (Spain) (Boselli & Gavazzi 
2002), at San Pedro Martir (Mexico) (Gavazzi ct 
al. 2002) and at the INT (Boselli et al. 2002). The 
estimated error on the Ha+[NII] flux is <~ 20%. 
Another important ingredient in this work, al- 
though it does not enter into the SEDs determina- 
tion, is the estimate of a galaxy current neutral hy- 
drogen content. HI data are taken from Hoffman 
et al. (1996, and references therein). HI fluxes 
are transformed into neutral hydrogen masses with 
an uncertainty of ~ 10%. From the hydrogen 
mass the HI deficiency parameter (Dcffji) is com- 
puted according to Giovanelli & Haynes (1985): 
Def ffl =< logM ffl (T obs ,D°£ t s ) > -logMgf where 
the observed HI mass is compared with the value 
expected from an isolated (i.e. free from external 
influences) galaxy having the same morphological 



type T obs and optical linear diameter D° b ^ (for de- 
tails, see Haynes & Giovanelli 1984). 

4. The galaxy SEDs 

Out of the 125 obtained spectra, 119 (2 Seyfert 
or Liners and 4 unclassified galaxies arc not in- 
cluded) are given in Fig. 3, grouped in 10 bins of 
Hubble type. These template spectra of normal 
cluster galaxies allow us to trace the dependence 
of the mean spectral properties along the Hubble 
sequence. Red continua characterize early type 
galaxies up to Sab-Sb (included), while the con- 
tinua become progressively bluer for later types. 
Except for dEs, which have weaker absorption 
lines, also the line properties of galaxies up to 
Sab-Sb appear indistinguishable, including the 
4000 A break. Emission lines are absent among 
galaxies earlier than Sa. Sa and Sab-Sb show no 
lines other than Ha and [Nil]. For later types the 
emission lines become progressively stronger, in- 
cluding [OIII], [Oil] and the other Balmer lines. 
The accurate analysis of the line properties is post- 
poned to Paper II of this series. 
The method we use to constrain a galaxy star 
formation history consists of running B&C mod- 
els with a broad grid of model parameters, and 
fitting the synthesized spectra to the observed 
ones to derive a set of best-fit parameters. How- 
ever, given the limited wavelength coverage of the 
optical spectroscopy, a similar approach would 
not sufficiently narrow down the parameter space. 
Much stronger constraints can be derived if the 
spectroscopic data are combined with photometry 
taken over the broadest possible wavelength base- 
line. We thus combine our 125 spectra ("spectro- 
photometric sample") with UV, optical (UBV) 
and infrared (JHK) photometric data. For this 
purpose magnitudes are converted into fluxes, 
then normalized at 5500 A, as for the spectra. 
Before the combination, both spectroscopic and 
photometric data are corrected for reddening ac- 
cording to the prescriptions of Appendix B. In 
this way we obtain SEDs that cover the domain 
2000-22000 A. For 220 objects in the "photomet- 
ric sample" that are not part of the "spectro- 
photometric" one the B&C models are fitted to 
the photometric data alone, as in Bell & de Jong 
(2000). 

In order to extend our SEDs below 2000 A we 
include the estimate of the ionizing flux below 
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Fig. 3. — Template optical spectra in 10 bins of Hubble type. The number of averaged spectra in each bin 
is given in parenthesis. Note that the flux scale differs between the left and the right panels. 
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912 A, as obtained from the Ha flux. This 
is based on the assumption that the flux in the 
Balmer line Ha is due to recombination of hydro- 
gen atoms within the HII regions excited by < 912 
A photons from the stellar radiation field. The 
details of this calculation are given in Appendix 
A. 

An example of such an extended composite SED 
(for the galaxy VCC 1205) is shown in Fig.4. The 
dots with error bars represent the available broad- 
band photometry with its uncertainty The opti- 
cal spectrum is represented by the thick contin- 
uous line. The dashed horizontal line in the far 
UV represents the < 912 A flux estimate derived 
from Ha. The thin continuous lines represent the 
Bruzual & Chariot models fitted to the data (see 
the next section). 

5. The Bruzual & Chariot models 

We use the 2001 version of Bruzual & Char- 
lot (1993) population synthesis models. These 
models provide the time (T) evolution of the syn- 
thesized spectra of galaxies characterized by an 
initial metallicity (Z), a Star Formation History 
(SFH) and an Initial Mass Function (IMF): 

F X (T)[Z,SFH,IMF] (1) 

Two star formation histories are considered: 

a) the exponential SFH: 

SFR(T exp , r exp ) = 1/r exp(-T/T) (2) 

which describes the exponential time evolution of 
a burst at T=Q. 

b) the " delayed-exponential" SFH or " a la Sandage" : 

SFR(T San , T S an) = T/t 2 exp(-T 2 /2r 2 ) (3) 

which mimics the SFH first proposed by Sandage 
(1986) (see Fig. 5). The temporal evolution of this 
family of functions is a delayed rise of the SFR up 
to a maximum, followed by an exponential de- 
crease. Both the delay time and the steepness of 
the decay are regulated by a single parameter r. 
We always assume a Salpeter IMF (a = 2.35 from 
0.1 to 100 M©; Salpeter 1955), and explore a pa- 
rameter grid in Z and r. Z is let free to vary 
from 1/50 to 2.5 Z© in 5 steps: 0.0004, 0.004, 
0.008, 0.02, 0.05. These values correspond to the 
initial metallicity grid computed by B&C and we 
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Z = 0.02 T E<= =13.0 t r „ p =20.00 
Z = 0.02 T IIp =3.0 t ew =1.80 
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Log X(A) 

Fig. 4.— The extended SED of the Sc galaxy VCC 
1205: the dots with error bars represent the avail- 
able broad-band photometry with its uncertainty, 
while the observed optical spectrum is represented 
by the thick continuous line. The dashed hori- 
zontal line in the far UV represents the < 912 
A flux estimate derived from Ha. The SED is 
fitted with B&C models with exponential SFH 
and solar metallicity (thin continuous and dashed 
lines). In spite of having chosen two very differ- 
ent combinations of ages and t (T=3 with r=1.8 
and T=13 with t=20 Gyr) the two solutions have 
similar probabilities. This example illustrates the 
degeneracy between T and r. 

have not tried to interpolate between them, r 
varies from 0.1 to 25 Gyr in 45 approximately log- 
arithmic steps. All models are computed at two 
epochs: T=5, 13 Gyr. 

Fig. 6 shows a sample of the synthetic galaxy spec- 
tra obtained assuming T=13 Gyr and Ts a n,Ts a n, 
our preferred choice for these two options (see next 
Section). In the top panel we explore the effects 
of changing the time-scale of the star formation 
process r while keeping the metallicity fixed at the 
solar value. A broad range of spectral properties 
is obtained, most noticeably when r < 5 Gyr. In 
the bottom panel we keep r fixed at 3.5 Gyr, and 
let Z vary over its full parameter range. In this 
case a much narrower spectral variety is obtained. 
The fitting of these synthetic spectra to the 
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Fig. 5. — The time evolution of a star formation 
event "a la Sandage" is plotted for 9 values of t 
from its origin (T=0) to 15 Gyr. This SFH is 
characterized by a steep rise followed by an expo- 
nential decay. It becomes progressively shallower 
and the position of the peak is delayed in time 
with increasing r. At the assumed age of 13 Gyr 
(vertical line) the SFH is rising for r > 13 Gyr. 

observed galaxy SED is carried out with a chi- 
squared minimization procedure that identifies the 
Z and t parameters for the most likely matching 
synthetic SED. Sometimes there is no clear min- 
imum in the distribution of chi-squared values, 
and the values flatten out until one reaches the 
edges of the parameter grid. In this case we as- 
sign lower or upper limits (as appropriate in each 
single case) to the values of the corresponding 
parameters. Details of the fitting procedures are 
presented in Appendix C. 

6. Results 

6.1. T vs. t degeneracy 

Fig. 4 shows the main limitation of our method 
as an age estimator: the SED of VCC 1205 is 
fitted with two B&C models with T exp ,T exp , so- 
lar metallicity and with very different ages (T=3, 
13 Gyr) and decay times (r=1.8, 20 Gyr), both 
models giving results consistent with the obser- 




Log A(A) 

Fig. 6. — Synthetic B&C galaxy spectra for a 
galaxy age T=13 Gyr and a star formation his- 
tory "a la Sandage" . In the top panel the models 
have fixed initial metallicity (Z=0.02) and cover a 
range of r values, while in the bottom panel the 
models have fixed t = 3.5 Gyr and cover a range 
of Z values. 

vations. Therefore, beside the well known age- 
metallicity degeneracy for a single stellar popu- 
lation (Worthey 1994), there exists another type 
of degeneracy that can significantly affect our at- 
tempts at deriving an accurate reconstruction of a 
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galaxy star formation history. In fact, because of 
the existing strong correlation between T and r, 
in spite of the broad wavelength coverage of our 
data, our method cannot disentangle the system 
age from r, but it provides only r/T. 
This is further illustrated in Fig. 7 where the ratios 
TSan/Tsan computed for two different ages (T=5 
and 13 Gyr) are plotted one against the other, 
showing that t/T is approximately constant while 
T varies by more than a factor of two. Due to this 
intrinsic weakness of the method we are forced to 
run B&C models with T fixed arbitrarily. We have 
chosen to use T=5 and 13 Gyr, but hereafter we 
show results only for T=13 Gyr, because the as- 
sumption of a common age of 5 Gyr for all galaxies 
appears to be at odds with the very existence of 
high redshift galaxies. 

6.2. The choice of SFH 

Having set T=13 Gyr we show in Fig. 8 a com- 
parison between the observed color-color (U — V 
vs. B — H) relation for the galaxies in our sam- 
ple and the grid of model colors obtained with the 
two adopted SFH (Exp, San). While both SFH 
reproduce without problems the reddest observed 
colors, the exponential SFH fails to reproduce the 
colors of the bluest galaxies, reaching at best a 
color U — V ~ 0.5 mag redder than observed, even 
assuming t—20 Gyr (quasi continuum SFH). The 
Sandagc SFH does a better job at reproducing the 
full range of observed colors within the explored 
range of r and metallicity, since it includes rising 
SFR at T=13 Gyr, thus allowing bluer colors. For 
this reason, and because it offers a more realistic 
time evolution of the galactic star formation his- 
tory, we adopt hereafter the Sandage SFH as the 
best single parameterization for the star forma- 
tion history of galaxies of all types. By allowing 
t to assume negative values the exponential law 
can also have a rising SFR, but this solution is 
less elegant than the Sandage law. Another way 
to create these blue colors is through starbursts 
and episodic star formation, which can have con- 
siderable effect especially on small galaxies. 

6.3. Fit to the template SEDs 

As discussed in the previous two sections, we 
select the Sandage SFH as the simplest and per- 
haps most realistic representation of the SFH of 
galaxies consistent with a fixed age of T=13 Gyr, 
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Fig. 7. — The value of the ratio Tsan/Tsan for 
the best fitting model with T—5 Gyr against that 
for the best fitting model with T=13 Gyr, fur- 
ther illustrating the degeneracy between T and 
r. Symbol shape represents Hubble type: ellip- 
tical galaxies are represented with squares, spirals 
with circles. Open and filled symbols differentiate 
galaxies in terms of HI deficiency (only for spirals): 
open symbols for galaxies with Dcfjji > 0.5 while 
filled symbols for galaxies with Dcfni < 0.5. Sym- 
bol size provides an indication of the best fitting 
model metallicity (size is decreasing with increas- 
ing metallicity). The typical uncertainty is given 
in the corner of the figure. 

and we let Z and r vary as free parameters in our 
fitting procedure. 

First we apply the fitting procedure to our tem- 
plates SEDs. These were obtained extending with 
photometric data the spectral templates shown in 
Fig. 3, and are therefore very robust determina- 
tions of the typical SED for each Hubble type. 
The results of the fit arc shown in Fig. 9. Good 
quality fits arc obtained with r increasing approx- 
imately monotonically along the Hubble sequence. 
It appears that the whole Hubble sequence can be 
modeled assuming increasingly delayed star for- 
mation histories of longer duration. Early type 
galaxies (dE, E, SO) have r < 3 Gyr, Sa-Sc have 
intermediate 4 < t < 6 Gyr, with uncertainties 
~ 1 Gyr. Galaxies in the 3 latest bins of Hubble 
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Fig. 8. — The observed (corrected for internal ex- 
tinction) U — V vs. B — H color-color distribution 
is compared with the corresponding color-color re- 
lation obtained from the B&C models (lines) for 
two star formation histories: the exponential one 
(top) and the one "a la Sandage" (bottom) (same 
symbols as in Fig. 7). 

type are consistent with constant or even rising 
SFHs. Im and BCD have lower limits r > 7 Gyr. 
The behavior of the metallicity along the Hub- 
ble sequence is more erratic. Approximately so- 



lar metallicities are found mostly in early-type ob- 
jects, except for dEs which show a sub-solar metal- 
licity. Slightly sub-solar metallicities are found 
across the entire range of Hubble types. Beside 
dEs, also Im+BCD galaxies show sub-solar metal- 
licities. Fig. 10 shows the iso-confidencc contours 
from the fitting procedure as a function of r and Z . 
Contours are plotted at 0.68, 0.85, 0.99 probabil- 
ity. The figure shows how r is better constrained 
than Z . Typical errors on r (see Sa, Sb, Sc) are 
0.9 Gyr. The uncertainty on the metallicity is typ- 
ically ~ 0.012, i.e. approximately one step of the 
metallicity grid. 

6.4. The SFH of elliptical galaxies as a 
function of luminosity 

The fit to template SEDs discussed in the pre- 
vious section has shown a significant variation of 
the fit parameters along the Hubble sequence. In 
this and in the following section we apply the fit- 
ting procedure to the individual SEDs, separately 
for Early and Late type galaxies. The fits were 
inspected one by one, but due to their large num- 
ber, they are not shown individually. 
We explore the dependence of the fit parameters 
on the NIR luminosity Lh , which traces the bulk 
of the stellar mass. This tight correlation between 
stellar mass and NIR luminosity has been dis- 
cussed by Gavazzi et al. (1996a) and Pierini et al. 
(2002) for luminous spiral galaxies, where HI mass 
does not dominate over the stellar mass, and by 
Zibctti et al. (2002) for elliptical galaxies. Thus 
the analysis presented in this and in the following 
section should shed some light on how the SFH of 
galaxies varies with their systemic mass. 
We reiterate that individual elliptical galaxies 
in the "photometric" and " spectro-photometric" 
samples are fitted with B&C models with fixed 
SFH (Sandage) and T=13 Gyr and the depen- 
dence of the derived parameters r and Z are ana- 
lyzed as a function of the system luminosity Lh- 
Fig. 11 shows the shallow dependence of r on the 
luminosity Lh for elliptical galaxies (the linear re- 
gression of this and of the following relations are 
given in Table 5). For a range of Lh spanning 4 
decades, r varies between 2 and 4 Gyr: i.e. these 
galaxies have experienced a "short" burst of star 
formation in their early history (see Fig. 5), with 
giant ellipticals being only twice "earlier" than 
dEs. 
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Fig. 9. — The (corrected for internal extinction) template SEDs , obtained averaging the photometry and the 
individual spectra grouped in 10 bins of Hubble type (same as Fig. 3) are fitted with B&C models assuming 
the Sandage SFH at T=13 Gyr, letting Z and r as free parameters (given in the labels together with the fit 
probability) . 
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Fig. 11. — The dependence of the fitted r on H 
band luminosity for elliptical galaxies (same sym- 
bols as in Fig. 7). Arrows mark upper and lower 
limits (see Appendix C). 



Fig. 10. — The confidence intervals as function of 
t and Z for the 10 Hubble type classes (same as 
Fig. 9) are given with contours drawn at 0.68, 
0.85, 0.99 probability. 



The symbol size in Fig. 11 increases with de- 
creasing metallicity of the best fitting model. This 
helps showing that, beside the shallow r vs. Lh 
relation, there is a more pronounced dependence 
of the model initial metallicity on the system lu- 
minosity, as shown in Fig. 12. At any given Lh, 
lower metallicity systems seem to have a lower r, 
due to the age-metallicity degeneracy. Summa- 
rizing, elliptical galaxies have experienced consis- 
tently "short" bursts of star formation "early" in 
their history. Their metallicity increases from sub- 
solar (1/5 Z©) for the dwarfs to solar for the gi- 
ants. 

6.5. The SFH of spiral galaxies as a func- 
tion of luminosity 

Figure 13 shows that the dependence of r on 
the system luminosity Lh is significantly steeper 
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Fig. 12. — The dependence of the fitted initial 
metallicity on H band luminosity for elliptical 
galaxies (same symbols as in Fig. 7). A small ran- 
dom number has been added to the metallicity in 
order to avoid superposition of points. 
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Fig. 13. — The dependence of the fitted r on H 
band luminosity for spiral galaxies (same symbols 
as in Fig. 7). Arrows mark upper and lower limits 
(see Appendix C). 

for spiral than for elliptical galaxies. Dwarf Irrs 
have t of 10 or more Gyr (increasing SFR), while 
the most massive spirals have r of approximately 
3 Gyr, consistent with elliptical galaxies. Once 
again in this figure the symbol size is inversely 
proportional to the metallicity of the best fitting 
model. While we observe a lower mean metallicity 
for any range of Lh than in the case of elliptical 
galaxies, spirals have a global dependence of the 
model metallicity on Lh not too different from the 
one discussed above for the ellipticals, as shown 
in Fig. 14. This dependence of metallicity on Lh 
confirms earlier findings by Bell & de Jong (2000). 
However, as far as t is concerned, these authors 
claim that the principal correlation is the one be- 
tween r and the central K band surface brightness, 
followed by a less significant r vs. Lk relation. 
Our data show a marginally better correlation of 
r vs. Lh than that of r vs. /i e . Notice however 
that we use the effective H band surface bright- 
ness /j, e instead of the central one used by Bell & 
de Jong (2000), and this might contribute to the 
discrepancy. We repeated the residual test of the 
Bell & de Jong and we find that the residuals of 
the r vs. fj, e relation correlate with L H slightly 
better than the residuals of the r vs. L H relation 



Fig. 14. — The dependence of the fitted metallic- 
ity on H band luminosity for spiral galaxies (same 
symbols as in Fig. 7). A small random number has 
been added to the metallicity in order to avoid 
superposition of points. 

correlate with /i e . Based on our analysis alone it is 
therefore impossible to discriminate whether /i e or 
L H is the principal parameter in regulating galaxy 
star formation history. 

6.6. The color-magnitude relation 

The color (B — H) vs. magnitude (Lh) relation 
for all galaxies, irrespective of their morphological 
type, is given in Fig. 15. Elliptical galaxies obey to 
a shallow color-magnitude relation which can be 
understood almost exclusively in terms of increas- 
ing metallicity with mass (the metallicity of dEs 
is significantly lower than that of giant ellipticals), 
with only a marginal spread in t in the explored 
luminosity range. The steeper (and non-linear) 
color-magnitude sequence of spirals/Im+BCD, in- 
stead, derives from the combination of increasing 
Z and decreasing t from dwarfs to giants. Notice 
that the colors of ellipticals are as red as those 
of spirals at the highest luminosities, suggesting a 
smooth transition from elliptical to spiral galaxies, 
with SO galaxies in between. 
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Fig. 15. — The (corrected) color-magnitude (B — 
H vs. Lh) relation for elliptical and spiral galaxies 
(same symbols as in Fig. 7). The dashed line is the 
fit for the elliptical galaxies. 



6.7. The environmental dependence of the 
SFH of spiral galaxies 

The results of the previous sections are prob- 
ably not unique to the Virgo cluster. They are 
consistent with Bell & de Jong (2000) analysis of 
non-cluster galaxies and it would not be surpris- 
ing to see similar luminosity dependencies among 
the isolated galaxies of Jansen et al. (2000). Are 
there features depending on the specific cluster en- 
vironment? Ideally one should investigate whether 
the r - luminosity relation varies as a function 
of some environmental parameter, like the local 
galaxy density or the projected angular distance 
from the cluster center. However Virgo has a com- 
plex structure, reflecting its non-virialized nature, 
which makes this approach unfeasible. Beside the 
main cluster A containing M87, a number of clouds 
are known to exist in this cluster: cluster B dom- 
inated by M49, clouds W, M, E and S, each with 
different velocity or distance, or both, as analyzed 
by Gavazzi ct al. (1999). As a result of the su- 
perposition of clouds on any line of sight, both the 
local galaxy density or the projected angular dis- 
tance from M87 provide little indication about the 
real environment surrounding any given galaxy in 



the Virgo cluster. 

We are therefore using a slightly indirect approach 
to the analysis of environmental effects, based on 
the Defni parameter, as defined by Giovanelli & 
Haynes (1985) (see also Sect. 3.3). Analyzing HI 
data of spiral galaxies in 9 rich clusters of galaxies, 
these authors concluded that HI deficiency is most 
likely due to ram-pressure stripping, a dynami- 
cal phenomenon occurring to galaxies in their fast 
motion through the cluster intergalactic medium 
(Gunn & Gott 1972). The extension of the sample, 
from the original 9 to 18 clusters, allowed Solanes 
at al. (2001) to conclude that the most HI defi- 
cient objects are those in radial orbits. The fact 
that asymmetries in the HI distribution (Dickey & 
Gavazzi 1991, Bravo-Alfaro ct al. 2000) are often 
associated with radio continuum trails (Gavazzi & 
Jaffc 1987, Gavazzi ct al. 1995) and in one case 
with an Ha trail (Gavazzi et al. 2001b) pointing in 
the same direction as the HI asymmetry, strongly 
suggests that ram-pressure stripping is the most 
plausible explanation for the pattern of HI defi- 
ciency in clusters. Thus the Defni parameter can 
be considered an indirect indicator of the environ- 
mental conditions experimented by a galaxy over 
the last few Gyr. A careful inspection of Fig. 13 
shows that galaxies suffering from significant HI 
deficiency (Defni >0.5) (empty circles) are found 
systematically below the regression line of the r 

- luminosity relation. The residual of the r - lu- 
minosity relation versus Defni is plotted in Fig. 16 
showing a significant trend. The deficient galaxies 
have SFH ~4 times shorter than the unperturbed 
galaxies. Summarizing: spiral galaxies have lumi- 
nosity dependent SFHs, and the residual of the r 

- luminosity relation is found to depend on the HI 
content, with the gas-poor galaxies showing trun- 
cated SFH compared with cluster galaxies with 
normal HI content. 

7. Discussion and conclusion 

We have combined spectro-photometry for 125 
galaxies in the Virgo cluster with broad-band pho- 
tometry from UV to NIR (including an estimate 
of the < 912 A flux from Ha). The resulting 
SEDs (corrected for internal extinction) were fit- 
ted with Bruzual & Chariot population synthesis 
models with a variety of model parameters (SFH, 
T, r, Z). The main results of the present work can 
be summarized as follows: 
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Fig. 16. — The residual of the relation between r 
and H band luminosity for spiral galaxies (Fig. 13) 
is given as a function of the HI deficiency param- 
eter (same symbols as in Fig. 7). 



1) the observations are consistent with the sim- 
plest (minimum number of assumptions) hypoth- 
esis that the all galaxies began their formation 13 
Gyr ago. This assumption implies that the star 
formation history of galaxies is allowed to increase 
with time, i.e adopting a SFH "a la Sandage" (see 
Fig.5). According to this family of SFHs the SFR 
increases with time, reaches a maximum and de- 
creases exponentially. A single parameter r deter- 
mines the delay of the SFR peak, thus the mean 
" age" of the stellar population, as well as the shape 
of its time dependence: the shorter r the more 
peaked the SFHs. 

Having made these assumptions the Hubble se- 
quence is obtained with increasing r with lateness, 
as proposed by Sandage (1986). 
2) r is a decreasing function of the H-band lu- 
minosity. This dependence is shallow for ellipti- 
cal galaxies: i.e. the value of r for the dwarfs 
(r ~ 4 Gyr) and that for the giants (r ~ 2 Gyr) 
do not differ by more than a factor of 2, both be- 
ing short compared with the assumed galaxy age 
(see Fig. 11). Late-type galaxies display a broader 
range of r values, from 25 Gyr (increasing SFR) 
for the low mass dwarf Im+BCD galaxies, down 
to 3 Gyr for the massive, early type spirals (see 



Fig.13). 

3) the stellar metallicity increases as a function 
of the H-band luminosity across the entire Hubble 
sequence (see Figs. 12 and 14). Z is found to span 
approximately a decade from low mass to massive 
systems. It should however be stressed that metal- 
licities are poorly constrained by the model. 
In short: the galaxy mass, as estimated from the 
H-band luminosity (or perhaps the surface bright- 
ness), seems to play a crucial role in determining 
most galaxy structural parameters (similar depen- 
dences of several other photometric and structural 
parameters on H-band luminosity are discussed in 
Scodeggio et al. 2002). 

Trying to reconcile the above results with mod- 
els of galaxy formation-evolution is beyond the 
scope of the present investigation. We refer the 
reader to Elmegreen (2001) for a comparison of the 
predictions of the hierarchical and of the mono- 
lithic models with the observations, unfortunately 
not explicitly focused on the dependence of the 
model predictions on the galaxy mass, an issue 
that we would like to address in the remainder of 
the present discussion, at least qualitatively. 
The hierarchical scenario contains an implicit ref- 
erence to the mass. The mass in this model is a 
growing function of time (i.e. of the number of 
mergings) . It is natural that the most massive ob- 
jects (either giant ellipticals or bulge-dominated 
spirals), being the product of a small number of 
major merging events which were likely to take 
place at z > 2 (Lacey & Cole 1993), can be mod- 
eled in our formalism with r ~ 3 Gyr. Low mass 
Im+BCD with significant current SFR, can be 
reconciled with the hierarchical scenario because 
dwarf galaxies either presently in formation or cur- 
rently undergoing minor mergings are allowed by 
the model. Their SFH, which we model for sim- 
plicity with shallow, continuous functions, may in 
fact consist of a series of minor star formation 
bursts, associated with minor merging events (see 
Kauffmann, Chariot & Balogh 2001). 
Our findings can in principle be reconciled also 
with the monolithic scenario if the mechanism of 
proto-galactic collapse is assumed to be not scale- 
free: i.e. that it depends on the total mass of the 
proto-galactic cloud, beside the dependence on its 
angular momentum. Angular momentum would 
discriminate between the formation of ellipticals 
and spirals with identical mass (Sandage 1986), 
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whereas the system mass would further regulate 
the efficiency of collapse: massive proto-galaxies 
would undergo an early efficient collapse, while 
small primordial fluctuations would take a longer 
time to form galaxies. As discussed by Boselli et 
al. (2001), not only the star formation history, 
but also the time dependence of the gas consump- 
tion can be modeled within the monolithic sce- 
nario with the above assumption. 
The additional finding of the present work is that 
the star formation time-scale r of spiral galaxies in 
the Virgo cluster depends on their present hydro- 
gen content. Highly HI deficient spirals in Virgo 
have typical time-scales (r) up to 4 times smaller 
than their HI healthy counterparts. If, as proposed 
by Giovanelli & Haynes (1985) and confirmed by 
Solanes et al. (2001), HI deficiency occurs to clus- 
ter spirals due to ram-pressure stripping (Gunn 
& Gott 1972), the exhaustion of the HI reser- 
voir might have produced an earlier truncation of 
the star formation activity in these galaxies. This 
mechanism is invoked by Gavazzi et al. (2002) to 
explain the significantly quenched current star for- 
mation rate found in HI deficient Virgo spirals in- 
cluded in their extensive Ha imaging survey. One 
might argue that the argument could be reversed: 
these galaxies appear HI deficient because a higher 
gas consumption occurred during an intense ear- 
lier star formation burst. However this reversed 
argument is in contradiction with the observations 
because it would imply that HI deficient galaxies 
would exist also outside clusters, something that 
is not observed. 

Among the highly HI deficient objects with trun- 
cated SFH there are few bright Sa in our sam- 
ple (VCC958 = NGC4419, VCC984 = NGC4425, 
VCC1047 = NGC4440, VCC1158 = NGC4461, 
VCC1412 = NGC4503) with r as short as 3 Gyr. 
This look-back time corresponds to z ~ 3, suggest- 
ing that at this early cosmological epoch Virgo was 
a fully developed cluster with a dense IGM. This 
is uncomfortably close to the maximum time of 
cluster formation allowed by the hierarchical the- 
ory (z ~ 2). 
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A. Estimate of the < 912 A flux from Ha + [Nil]. 

The stellar radiation field with A < 912 A ionizes the gas which re-emits, via recombination lines. If the 
gas is optically thick in the Lyman continuum, the number of photons in a specific recombination line is 
directly proportional to the number of star photons in the Lyman continuum. In the case of H/3 this number 
is given by equation (5.23) in Osterbrock (1989). For Ha we have: 



oc 

/ 



^du = L Ha -C (Al) 



where: 

Assuming T=10000K and the Osterbrock case B: 
a e JjH(H°,T) = 3 x l(T 14 (cm 3 sec" 1 ) 
a B (H°,T) = 2.59 x l(r 13 (cm 3 sec" 1 ) 
and ^ = 2.87 

From the H Q luminosity it is thus possible to recover the number of ionizing photons, which can be compared 
with the similar quantity derived from the integral on the model spectrum. 

For galaxies in the " spectro-photometric" sample Ln a is measured. For galaxies in the "photometric" sample 
Lu a can be in most cases derived from imaging observations. These however must first be corrected for the 
contamination of [Nil]. From our spectra we calibrate an empirical relation between Ha/ [Nil] and Lh and 
apply it to estimate [Nil] for the remaining galaxies of known Lh- 

The derived Ln a must be further corrected for extinction. We assume 0.8 and 0.6 mag extinction at Ha for 
galaxies earlier than Scd and for Sd-Im-BDC respectively (consistently with Boselli et al. 2001). Moreover 
we assume that 100% of the UV photons contribute to excite the recombination lines (Mezger 1978), with 
en escaping fraction of 0. A conservative estimate of the uncertainty on the derived < 912 A flux is 1 mag. 



B. The extinction correction 

The issue of determining the internal extinction in galaxies is a very controversial (from the extreme 
optically transparent assumption of Holmberg, 1958 to the optically thick case of Valentijn 1990) and difficult 
one, due to the complexity of the mutual geometry of obscuring dust and stars (Chariot & Fall 2000). 
Moreover stars of different ages have different scale heights, making the determination of the wavelength 
dependence of extinction a very uncertain task. The scale height ratio dust / stars was found to vary with A 
from 0.2 in H to 0.9 in B by Boselli & Gavazzi (1994), whereas Xilouris et al. 1999 found it approximately 
constant. 

Solving this issue is beyond the scope of the present investigation. However we cannot avoid correcting our 

extended SEDs for internal extinction using some, even simplistic recipe. 

The extinction corrected intensity at any A is related to the observed intensity by: 

LogI xc = Logh + 0.4 • A\ (Bl) 

where 

Ax = f(r x ) (B2) 

t(A) is the optical thickness at A and / depends on the relative geometry of dust and stars. 

We tested several approaches, whose results are shown in Fig. 17, carrying the SED of the moderately 

inclined Sc galaxy VCC 2058, before and after the correction for internal extinction. 

1) We tried the Calzetti (2001, C01 hereafter) empirical absorption law derived for star-burst galaxies from 
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the Balmer decrement (top-right panel of Fig. 17). As pointed out by Buat et al. (2002) the Calzetti law 
can hardly apply to normal galaxies, as it overestimates the absorption in the UV, due to the particular 
geometry of star-burst galaxies. 

2) We derived A uv = g(FIR/UV) from the FIR/UV ratio, as proposed by Buat et al. (1999, 2002) and 
Witt & Gordon (2000; W&G00 hereafter). This choice is supported by the fact that FIR and UV fluxes are 
in first approximation independent from the SFH (the two are produced by similar stellar populations) and 
from the assumed dust-star geometry. From Auv we calculate tjjv for 3 geometries, than we derive 

rx = I'(Ax) (B3) 

tx = r uv ■ kx/k uv (B4) 

using the galactic extinction law k{\) (Savage & Mathis 1979), and we compute the complete set of Ax using 
(B2). 

The assumed geometries are: 

2a) the sandwich model of Disney et al. (1989, D89 hereafter) (bottom-right panel of Fig. 17). We used a A 
dependent dust to stars scale ratio: 0.86(B), 0.78(V), 0.18(H) (intermediate between the optically thick and 
thin case of Bosclli & Gavazzi 1994). 
2b) the slab model of D89 (bottom- left panel of Fig. 17). 

2c) the dusty-spherical slab model (with clumpy dust distribution) of W&G00, which includes the contribu- 
tion from the scattered light (top- left panel of Fig. 17). 

As concluded earlier, we rule out the Calzetti law because it overestimates Auv in normal galaxies. Sim- 
ilarly we excluded the W&G00 model because the spherical approximation is hardly a realistic one. Since 
the sandwich and slab models of D89 produce very consistent results we finally adopted for simplicity the 
slab model. 

In this model the explicit functional dependence of (B2) gives the face-on: 

Ax = 2.5 • Log l - ^P(-T, ■ secW 
tx ■ sec(i) 

where sec(i) = a/b. According to Buat et al. (1999) Ajjv can be quantified as: 

A uv (mag) = 0.466 + Log(FIR/UV) + 0.433 x (Log(FIR/UV)) 2 (B6) 

where FIR = 1.26 x (2.58 x 10 12 x F 60 + 10 12 x F WQ ) x 10~ 26 [Wra' 2 ] 

F 60 and F wo are the IRAS FIR fluxes (in Jy) 

and UV = 10- 3 x 2000 * io(^m. 9 +2i.i75)/-2.5 [Wm~ 2 ] 

(B4) can be obtained from a polynomial fitting inverting (B6): 

tuv = (l/sec(i)) x (0.0259 + 1.2002 x A uv + 1.5543 x A 2 UV - 0.7409 x A uv + 0.2246 x A\, v ) (B7) 

Thus Ax are finally obtained. FIR/UV is available for 85 objects. If FIR or UV measurements are unavailable 
we assume Auv — 1-28; 0.85; 0.68 mag for Sa-Sbc; Sc-Scd; Sd-Im-BCD galaxies respectively, as determined 
for galaxies in these three classes of Hubble type when FIR and UV measurements are available. 
After the SEDs were corrected with the above method, we checked empirically that they do not contain a 
residual dependence on the galaxy inclination. The corrected SEDs of 32 Sc galaxies, binned in 4 intervals 
of inclination, and their fit parameters were found very consistent one another. 

C. Fit of B&C models to the data 

The logarithmic fit of B&C models to the data is obtained minimizing the \ 2 (reduced chi-squared): 



* 2 = l X 



-^i,obs ^i,'mod 



d . 
i=i 



(Cl) 
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where d is the number of degrees of freedom, d = n — v, n represents the number of bands for which i*i )0 j, s =/= 0, 
v is the number of constraints, i.e. the number of fitted parameters (e.g. Z and r). Fi,obs and Fi, mo d are the 
observed and model fluxes respectively, and <r, is the error associated with each photometric "band" (UV, 
U, B V, J, H, K, Ha, spectrum). 

logF hobs = -0.4 • Ma 9l + C t - A(f t ) (erg cm" 2 sec" 1 A" 1 ) (C2) 

where 

d = -2.5 • (log(J F Vega (X)R t (X)d\) 
are tabulated in Table 6 and 
A(fi) = -2.5 -logiJ Ri(X)dX) 

represent the area of the filter transmission (taken from B&C). 

For the B, V, H bands we assume an uncertainty of 0.15 mag, 0.2 for U, J and K, 0.4 mag for UV and a 
large uncertainty of 1 mag for the < 912 A flux derived from Ha. 

The spectra were divided in 6 bins: 3800-4000; 4000-4200; 4200-5000; 5000-5800; 5800-6600; 6600-7000 A. 
The first 2 narrow bins bracket the 4000 A break. Each spectral bin contributes to the fit as one "photometric 
point" . The error on each point was determined by combining the rms measured in line-free regions of the 
individual spectral bins with 0.10 mag uncertainty in the photometric calibration. 
From x 2 we compute a fit probability defined as: 

p( * 2 ^ = ^rpj k/~ le ~ x2 ' 2dx - (C3) 

where V is the incomplete-gamma function. P(x 2 > xlbs) is the probability of obtaining a x 2 equal to the 
observed value {Xobs)- Only fits with probability > 5% are considered. Only 3 galaxies in the photometric 
sample and none in the spectro-photometric sample have P < 5%, while 85 % of the objects have P > 70%. 
For any given metallicity the probability of the fit is traced as a function of r. For most cases the distribution 
reaches a maximum at a given t and Z which are then selected as best parameters (see Fig. 18 center). In 
some cases the probability never reaches a maximum within the parameter space, but tends asymptotically 
to a maximum either for t — > or t — > oo (see Fig. 18 left and right). Because neither of these values is 
meaningful in these cases we adopt t, Z corresponding to a threshold probability of 0.90 of the maximum. 
Two notes of caution. 

1) In spite of our effort for homogenizing the apertures at which data used in the present work were taken, 
total UV magnitudes, integrated spectra and broad band photometry have still slight aperture differences 
(see Section 3). The residual aperture effects on the SED shapes should however be within the quoted 
photometric uncertainties, thus not affecting significantly the determination of the fit parameters. 

2) The incompleteness in the various observations bands increases toward lower luminosities (see Section 2). 
This might bias the luminosity dependence of the fit parameters. We have checked this effect by excluding 
one by one from the fit the UV, the U and B data and the spectra. We conclude that none of these exclusions 
produces significant changes. In particular excluding the U, B and spectral data does nothing but adding 
some noise in the correlations. The exclusion of UV data has an appreciable effect on the SEDs of some 
early- type spirals which become "bluer", thus with longer r. 
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Filter 


C 


UV 


15.44 


U 


14.11 


B 


13.06 


V 


13.77 


J 


15.85 


H 


16.11 


K 


16.62 



Table 6: Conversion coefficients from magnitude to flux for the considered filters. The coefficient for the UV filter corresponds 
to 2000 A since UV1650 magnitudes are converted to 2000 Ausing an average correction of 0.2 mag (see Section 3.2). 




Fig. 17. — The SED of VCC2058 is shown uncorrected (light curve) and corrected according to the dusty- 
spherical slab model of Witt & Gordon (2000) (top-left panel); Calzetti (2001) (top-right panel); the slab 
model of Disney ct al. (1989) (bottom-left panel) and the sandwich model of Disney et al. (1989) (bottom- 
right panel). 
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Fig. 18. — The fit probability as a function of r for the various metallicities. An example of P(r) reaching 
a maximum for t — > is given by VCC 758 (left). An example of P(r) reaching a maximum within the 
parameter space is given by VCC 66 (center). In this case r, Z max are univocally determined from P ma x(T, Z). 
An example of P(r) reaching a maximum for r — > oo is given by VCC 159 (right). In the first and third case 
t, Z max are determined using the thresholded probability of 0.90 of their maximum. 
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